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■ ABSTRACT 

Context. The PICsIT detector onboard the INTEGRAL satellite was designed to provide information about emission in the soft y-ray 
band for many bright sources. Due to strong and variable instrumental background, only 4 objects have been detected so far using 
standard software. 

Aims. The moderate sensitivity of PICsIT can be compensated for in the case of many objects by adopting a long exposure time, thanks 
to integral's large field of view and an observing strategy focused on the Galactic plane. With angular resolution far higher than 

■ that of all other instruments operating in a similar energy band, PICsIT is suitable for fields too crowded or too significantly affected 
by Galactic diffuse emission. Therefore, it is desirable to improve the spectral extraction software to both obtain more reliable results 
and enlarge the number of objects that can be studied. 

Methods. The new PICsIT spectral extraction method is based on three elements: careful modelling of the background, an energy- 
dependent pixel-illumination function, and the computation of the probability density of the source count rate. Background maps of 
^ ' the detector plane are prepared for short periods and relatively narrow energy bands to insure that the background dependence on 

' time and energy is modelled well. The most important element of the new spectral extraction method is the proper treatment of the 

Poisson-distributed data, developed within a Bayesian framework. 
^S) ■ Results. The new method was tested extensively on both a large true data set and simulated data. Results assumed in simulations 

^ ■ were reproduced perfectly, without any bias and with high precision. Count rates measured for Crab were far more stable than those 

' obtained with the standard software. For weaker sources, the new method produced spectra of far higher quality and allows us to 

I detect at least 8 additional objects. Comparison with other INTEGRAL instruments demonstrated that PICsIT is well calibrated and 

, provides valuable information about the continuum emission in the 250 keV - 1 MeV band, detectable currently only by INTEGRAL. 

^ , Key words. Gamma rays: observations - Instrumentation: detectors - Methods: data analysis - Methods: statistical 
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99 1 . Introduction 13', and its nominal energy range was 35-1300 keV. The OSSE 

T\ telescope has been the most sensitive soft y-ray instrument sent 

^ The low-energy y-ray band (100 keV - 10 MeV) is extremely -^^^ ^^^-^ ^^^^ -^^ j^^g^ ^-^^ ^j^^ j^igj^ sensitivity of 

■ difficult to observe mainly because of the intense nuclear-decay Q^^g p^^^^i^^ ^j^gj^^j^ ^^^^^^^ limitations related to 



■ 

X 



emission induced by cosmic rays inside the detector and sur- -^^ ^^^-^ ^^^-^-^^^ j ^^j^ ^-^^ qSSj, ^^^^ 

H . rounding materials. As the flux of the observed objects decline ^^^^^^ ^^^^j^^ ^j^^^^ ^^-^^^^ ^1^^ q^^^^^-^ ^^^^^^ 

. ^ , significantly witii mcreasing energy from X-rays to y-rays, the distinguish point-source emission from Galactic diffuse emis- 

strong increase in internal background of the detector makes ^-^^ g^^^g ^jj.^^ ^^^-^^^^ ^^^^ ^^-^^ detecting 

the ratio of observed signal to background very small. This ra- q^.^^ monitoring the activity of bright er y-ray source s. The 

tio usually does not exceed 0.01 even for sources as bright as ^^^-^^^ pi^S dFrontera et alJll997h 

the Crab nebula, despite the passive and active shielding ap- ^^^^^^ BeppoSAX satellite was 15-300 keV but its sen- 

phed. Moreover, the fact that soft y-rays cannot be focused as ^-^^^ ^^^^^ insufficient to provide more precise 

well as soft X-rays reduces the angular resolution in this en- ^^j^ 
ergy range. Therefore, observations of crowded Galactic fields 

with y-ray telescopes are often affected by contamination from There are four missions presently in operation with soft y- 
nearby sources or from Galactic diffuse emission. ray instruments, namely RXTE (HEXTE), IN TEGRAL, Swift 
There have been many y-ray observatories, but due to their (BAT ), and Suzak u (HXD/GSO). The HEXTE dRothschild et al.1 
generally limited life times and sensitivities only a few have pro- 1998') and BAT OGehrels eta l.l l2004 energy ranges do not ex- 
vided a significant amount of data. More detailed studies of the ceed «150 keV, and these instruments have therefore limited ca- 
soft y-ray spectra became possible in the l ast decade of the 20th pabilities of studying the high energy emission. GSO onboai'd 
cenhiry, thanks to the SIGMA detector (iPaul etal.llT99l on- Suzaku was declai-ed to be the detector with the lowest inter- 
boar d the GRANAT satellite, and the OSSE and BATSE detec- nal background in the 40-600 ke V band, due to a mo re ad- 
tors dJohnson et al .''l993) onboard the CGRO satellite. SIGMA vanced active-shielding technique dTakahashi e t al.l l2007l) . This 
was a coded-mask instrument with a good angular resolution of low background should correspond to high sensitivity, but 3 
years after Suzaku's launch almost no results above 200 keV 
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keV dMakishima et al.|[2008l) . Therefore, taking into account the 
hmited angular resolution of GSO above 100 keV 4.5°), cur- 
rent interest in the soft y-ray research is focused on the data ac- 
quired by INTEGRAL detectors^ 

The INTEGRAL satellite jWinkler et a l."2003') has onboard 
three coded-mask y-ray detectors. ISGRI (Lebrun et al. 2003) is 
the u pper detector of INTEGRAL'S imager IBIS (Uberti ni et al] 
12003 ). which operates nominally in the 13-1000 keV range. 
Owing to its relatively high sensitivity and angular resolution 
of 12', ISGRI remains the INTEGRAL instrument of choice for 
studies of point sources, particularly Galactic ones. The lower 
layer of the IBIS imager, the PICsIT detector (Laba nti et alj 
l2003h. is described in the next section. The spectrometer SPI 
dVedrenne et"ani2003l) is a high spectral resolution y-ray tele- 
scope operating in the energy range 20 keV - 8 MeV. Despite 
its lower sensitivity and limited angular resolution {a 2.5°), SPI 
spectra extracted for brighter sources complement the ISGRI 
spectra well, especially above 150 keV where the ISGRI sen- 
sitivity declines rapidly. Due to its reliable, ground-based effi- 
ciency calibration, SPI serves also as a reference instrument for 
the INTEGRAL cross-calibration tests. Observations of crowded 
fields, however, are difficult, a s for OSSE, and special ca re is re- 
quired in analysing such data (i Roques & Jourdainll2005l) . 

Notwithstanding their unique capabilities, SPI and ISGRI 
provide rather limited information about the continuum emission 
from point sources above «200 keV. In the case of SPI, besides 
its limited angular resolution, there are some difficulties in mod- 
elling the instrumental background at higher energy, where the 
number of photons is low. However, by applying a more sophis- 
ticated analysis to SPI data, it w as possible to detect 20 objects 
above 200 keV (B ouchet et all 12008) . The sensitivity of ISGRI 
declines dramatically above 150 ke V, and the efficiency c alibra- 
tion remains uncertain in that range (I jourdain et al.ll2008l) . 

Initially, the INTEGRAL imager sensitivity was pl anned to 
excee d that of OSSE by about an order of magnitude dWinkleij 
|1994|) . Unfortunately, due to a reduction in mission funds, this 
plan became impossible, when the foreseen three layers of the 
high energy detector were reduced to two layers of smaller vol- 
ume and a limited telemetry was ascribed to them. The effect is 
that the PICsIT detector, being still the largest soft y-ray detec- 
tor on orbit, has the thickness of only 3 cm compared to almost 
18 cm of OSSE. The sensitivity of PICsIT is reduced further be- 
cause the pixellated structure and data-acquisition logic produce 
an additional loss in the fraction of the Compton-scattered pho- 
tons. On the other hand, due to the elongated orbit, both the Earth 
occultations and passages through the radiation belts affect only 
a small part of the INTEGRAL revolution period. This and the 
observing strategy focused on the Galactic centre and Galactic 
plane ensures that many objects are observed by PICsIT with a 
long exposure time, compensating to some extent the moderate 
sensitivity of the instrument. 

Shortly after launch, it appeared that the PICsIT background 
level was about two times lower than the pre-launch conservative 
estimate (,Di Cocco et al. 2003.) . Only the energy band below 300 
keV showed a higher background, with an excess originating in 
track events induced by high-energy cosmic rays ("Segre to et alj 
l2003h . Despite the lower background, the sensitivity achieved af- 
ter the first year of detector operation was low, a factor between 
5 and 10 lower than the statistical limit. The situation was im- 
proved considerably when the long-exposure b ackground map s 
were included in the OSA 4.0 software release (lFoschinill2004D . 
allowing a sensitivity comparable to that of SPI to be achieved. 

Nonetheless, almost six years after the INTEGRAL launch, 
there remain only three objects for which results with PICsIT 



have been reporl:ed: C rab dDi Cocco et all l2003h . Cyg X-1 
dCadoUe Bel et al."2006'), and XTE J 1550-564 during its 2003 
outburst (Foschini 2005). The PICsIT spectra presented did 
not provide significantly more information than the SPI spec- 
tra. More spectacular PICsIT results were instead obtai ned for 
the o bservations of gamma-ray bursts (GRBs) (Malag uti et al.l 
l2003 a). in particular those appearing outside the INTEGRAL 
field of view (Marcinkowski et al. 2006), because these events 
are detectable at high energy with the use of IBIS Compton- 
mode data. 

Since the possibilities of improving the spectral extraction 
performance for PICsIT using the standard OSA software were 
limited, a completely novel approach was developed for test- 
ing purposes. Preliminary results demonstrated that this new ap- 
proach produced very good results, by for example allowing us 
to detect several objects not detected by the standard procedure. 
This paper presents all the elements of the new method, which is 
now fully developed. Since we decided after discussions with the 
IBIS Team that the new method should not be implemented in 
the standard software for INTEGRAL data analysis, we provided 
a detailed presentation here to allow a potential user to repro- 
duce the results. In Sects. |2]3 basic ingredients of the spectral 
extraction technique are described. Section|5]presents the results 
of testing the PICsIT detection limits. In Sect. |6l our main re- 
sults, PICsIT spectra, are widely presented and discussed. After 
a summary given in Sect. |2l an extended Appendix A presents 
the results of tests completed for different methods that can be 
used to extract the source count rates from y-ray instruments. 

2. Instrument model 

A detailed description of the PICsIT detector an d IBIS in genera l 
can be found in the IBIS Observer's Manual ( lKuuIkersll2006h 
and in the references listed there. We present here only basic 
information related to the pixel-illumination modelling. 

PICsIT is the Pixellated Caesium Iodide Telescope, which 
operates in the nominal energy range of 175 keV - 10 MeV and 
has an angular resolution of 12'. The fully coded field of view is 
a rectangle of 9°x 9° and the 50% partially coded field of view 
has a size of 19°x 19°. The detector consists of 4096 Csl crys- 
tals, organized in 8 Modular Detection Units. Each pixel has the 
surface area of 8.55x8.55 mm^ and a 30 mm thickness; pixels in 
a module are separated by 0.55 mm. The total sensitive area is 
2994 cm^. The IBIS mask is located 3283 mm above the PICsIT 
surface. The mask consists of 9025 open or closed cells, each 
of surface area 11.2x11.2 mm^, the total size of the mask be- 
ing 1064x1064 mm^. Closed cells consist of 16 mm thickness 
of tungsten. The ISGRI detector made of 16384 CdTe pixels 
(4x4x2 mm-') is mounted 9 cm above the PICsIT surface. The 
IBIS collimator operates as a passive lateral shield and consists 
of the main IBIS tube, a hopper giving an additional limit to 
solid angle observable by detectors, and a side-mask shielding 
between the mask and tube walls. Both ISGRI and PICsIT de- 
tectors are surrounded with lateral and bottom BGO crystals that 
form an anti-coincidence shield, reducing the instrumental back- 
ground by about a factor of 2. 

PICsIT operates in three modes: photon-by-photon, spectral- 
imaging, and spectral-timing. Due to the limited INTEGRAL 
telemetry, photon-by-photon mode is in practice used only for 
tests. Therefore, almost all PICsIT data are stored as spectral- 
imaging histograms and spectral-timing tables. The spectral- 
imaging is the basic mode since it provides full position infor- 
mation (count rate for each pixel) and 256 channels of energy 
information. Histograms correspond to a single science window 
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(a single spacecraft pointing), thus the time resolution of the 
spectral-imaging data is usually in the range of 30-120 minutes. 
Spectral-timing data have a time resolution fixed between 1 and 
500 ms but the count rates are integrated over all detector pixels 
and stored in only 8 energy channels. Therefore, this mode can 
be used only for extremely bright sources, typically GRBs, that 
have the signal strength that exceeds the background emission in 
short time intervals. 

The spectral-imaging histograms are made for two types of 
events: single and multiple. A single event detection is one in 
which there is only one photon in a single PICsIT pixel. Multiple 
events are defined to be detections of 2 or 3 photons in a single 
submodule (half of the module, with pixels coupled electron- 
ically); the event position is then given by the position of the 
pixel with the highest energy deposit and the energy of the event 
equals the summed energy of all photons registered within the 
event time. The fraction of multiple events varies from x 10% 
in the 300^00 keV band to ^ 50% in the 1500-2200 keV band. 
The logic of selecting single and multiple events based on one 
submodule detection introduces some non-uniformity because 
the multiple events in the border pixels are often registered as 
single events. This paper describes the spectral extraction ap- 
plied to single events only, although the same procedure can be 
applied to multiple events. 

The standard analysis software for all INTEGRAL data is 
the Off-line Science Analysis (OS A) package, distributed b y the 
INTEGRAL Science Data Centre jCourvoisier et al.ll2003h . The 
OSA software package presently allows for a limited PICsIT 
data analysis: image analysis for spectral-imaging data, and de- 
tector light curve and spectrum analysis for spectral-timing data. 
Taking into account that each source detected by PICsIT should 
be observable at lower energies by ISGRI, PICsIT alone is not 
expected to detect new sources. Therefore, the spectral extrac- 
tion based on the Pixel Illumination Function (PIE) computed 
for the catalog source position appears to be more appropriate for 
quantitative PICsIT data analysis than the image deconvolution. 
The PIE-based method directly describes the photon absorption 
that on the detector surface creates a shadow of the mask illumi- 
nated by the source observed at a given angle. In contrast, image 
deconvolution is an indirect method, which usually applies some 
type of inversion or cross-correlation that can lead to decreased 
efficiency. 

There has been a spectral extraction component in the OSA 
software sinc e version 5.0 but the PIE model included is not yet 
reliable (Eos cliinil[2"007al) . In consequence, all PICsIT spectra 
made with OSA are based on count rates determined from the 
mosaic images. Eor the purpose of the research presented here, a 
completely new software was developed for PIE modelling. The 
code was an adaptation of the code used to model th e absorp- 
tion in the Nomex structure supporting the IBIS mask dLubiriskil 
|2007| ). The absorption of photons in the mask and the ISGRI de- 
tector layer was computed to be the mean fraction of source pho- 
tons arriving at each pixel at given polar and azimuthal angles. 
Attenuation coefficients for tungsten (mask) and CdT e (ISGRI) 
are taken from the tables of iHubbell & Seltzej ( Il996h . The PIE 
map is calculated for 22 energy bins, with an effective energy for 
each bin that is calculated by assuming that the spectral index of 
the source emission is 2.1. The probability that the photon ar- 
rives at each pixel of the detector is computed over a grid where 
photons are separated by 0.125 mm. The photon path length in- 
side a given mask or ISGRI pixel was determined by numerical 
integration and the vertical step was also set to the value 0.125 
mm, which was far smaller than the pixel size and thickness. 
Eigure [T] presents an example of ISGRI and PICsIT PIE mod- 



els calculated for the same science window and the same source. 
The Nomex structure supporting the IBIS mask absorbs between 
10% and 25% of the source photons at energies above 200 keV, 
depending on the source off-axis angle. To correct the source 
count rates for this effect, a set of high-energy off-axis correc- 
tion maps was prepared for PICsIT, in the same manner as those 
provided for ISGRI in the OSA package. 





Fig. 1. ISGRI (left, 50 keV) and PICsIT (right, 450 keV) PIE 
models for science window 007900330010 calculated for Cyg 
X-1 seen at 3.9° off-axis angle and 224° azimuthal angle. 
Azimuthal angle is measured anticlockwise from the right. The 
PICsIT PIE is slightly shifted towards the upper right because of 
larger distance between the mask and detector EuUy open and 
fully closed elements transparency is respectively 1 and 0.003 
for ISGRI and 0.89 and 0. 1 1 for PICsIT. 

High energy photons from the source can pass through the 
IBIS tube walls and illuminate the part of the detector that is 
not coded by the IBIS mask. Eor that reason the non-coded 
part of shadowgrams cannot be used to estimate the background 
level and the PIE is usually calculated only for the coded part. 
However, the ISGRI PIE model implemented in OSA takes into 
account different structures such as walls or hopper, providing 
a source illumination pattern for the entire detector plane. This 
PIE model can be adapted for PICsIT after applying shifts and 
changing the size of the pixels. Tube walls are quite transparent 
to photons in the PICsIT energy range, for example, for sources 
at an off-axis angle of 36° and with an azimuthal angle of 53° 
about 40% of 450 keV photons can pass through the wall, illu- 
minating the detector surface. Therefore, it appears desirable to 
use such data in attempting to increase the amount of informa- 
tion about a given source where the large off-axis angle obser- 
vations correspond to the exposure time, which can be several 
times longer than the time of fully coded observations. Eigure 
12] shows an example of the extended PIE model for the source 
observed at 36° off-axis. The extended PIE model for PICsIT 
spectral extraction was tested with the Crab data, and the results 
of these tests ai^e presented in Sect. 16.21 

3. Background model 

Given the fact that the PICsIT instrumental background rate ex- 
ceeds the source count rate by several orders of magnitude, a 
correct analysis of PICsIT data must be based on careful mod- 
elling of the background. The background model for a given 
observation must be prepared after detailed studies of general 
PICsIT background properties and investigation of the behaviour 
of the background during that observation. Positions of the y-ray 
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Fig. 2. PICsIT extended PIF model for science window 
0220009700 10 calculated for Crab seen at 36° off-axis angle and 
53° azimuthal angle. The mean wall transparency for 450 keV 
photons is 63% for the lower left part (red), 39% for the upper 
part (green) and 20% for the lower right part of the shadowgram 
(blue). 



lines are one of the basic characteristics of the background. Since 
these lines are emitted by isotopes produced by cosmic rays in 
the detector and the surrounding materials, a spatial distribution 
of the count rate over the detector surface is energy dependent. 
Another important property of the background is its time evo- 
lution. Long-term variability should be monitored to check, for 
example, the positions of the lines in the background spectra. 
Short-term variability must be studied for selection of data that 
can be used to prepare the background map for a given obser- 
vation. Solar flare periods or time intervals with a strange back- 
ground behaviour should also be excluded. This section presents 
the results of the PICsIT background studies. 

3.1. Spectra 

Figure |3] shows PICsIT background spectra from different peri- 
ods of the mission. The spectra were integrated over one space- 
craft revolution (hereafter Rev.) lasting almost 72 hours. The 
specti-um of Rev. 0541 (19-21 March 2007), taken exactly 4 
years after the Rev. 0052 spectrum (18-20 March 2003) illus- 
trates how the background emission evolves during the mission. 
The sharpest increase in the emission with time was observed 
at low energy, below 300 keV, where the Rev. 0541 count rate 
was about as twice high as the Rev. 0052 count rate. Above 300 
keV, the count rate after 4 years was about 50% higher than at 
the beginning of INTEGRAL operation. 

The PICsIT energy resolution is rat her coarse, varying be - 
tween 30% at 300 keV and 4% at 4 MeV (iMalaguti et al.l2003bh . 
In spite of this, the background spectra are not completely 
smooth, exhibiting strong peaks at 511 and 680 keV, weaker 
peaks at 1800 and 2750 keV, and a broad hump between 800 
and 1700 keV. The identification of the background characteris- 
tic features can be useful in properly adjusting the energy bins 
used in the data analysis and monitoring the energy calibration 
stability. Due to the limited energy resolution even the strongest 
lines positions were affected by contamination from weaker lines 
nearby, while weaker lines were hidden in the Compton con- 
tinuum associated with the stronger lines, which form broad 
humps everywhere below 1700 keV. Therefore, only several 
background lines could be identified with a high level of plausi- 
bility. These identifications were s upported by comparison with 
the list of SPI background lines jWeidenspointner et al.ll2003h 



Table 1. Strongest PICsIT background lines seen in the spectra 
collected during non-flaring Sun periods. The most probable line 
origin is listed first. T1/2 is the half-life decay time of a given iso- 
tope. Lines around 670 and 690 keV form the 680 keV complex 
(see text). 



Energy 
[KeV J 


Catalog energy 

[KeV J 


Isotope 


Tl/2 


510 


511.0 






670 


666.3 


126j 


13 d 




667.7 


"2CS 


7d 




666.3 




12 d 


690 


695.0 


12«Sb 


12 d 




685.7 


'2'Sb 


4d 




703.4 


205Bi 


15 d 


1780 


1779.0 


2'*A1 


2m 




1757.6 




36 h 




1771.4 


''Co 


77 d 


2750 


2754.0 




15 h 



that were clearly res olved, and the da ta about Csl proton activa- 
tion measurements dRuiz et al.lll994l) . The main lines identified 
in normal (non-flaring Sun periods) spectra are listed in Table 
[1] The two strongest lines, the 511 ke V (j6^ decay line) and the 
680 keV complex, are emitted predominantly by the iodine, cae- 
sium, and antimony isotopes produced in the detector material. 
Radioactive bismuth is produced in the BGO crystals of the IBIS 
anticoincidence system, whereas sodium and aluminum activity 
originates in the aluminum frame holding both IBIS detectors. 

In addition to the lines presented in Table [1] normal spec- 
tra also exhibit weak emission around 2230 keV, presumably 
from a 'H(n,7')^D reaction. The bumps observed between 800 
and 1700 keV are probably associated with emission from nu- 
merous isotopes produced mainly in the detector and its frame, 
for example ^^Fe, ^^Ni, ^'^Mg, and '^"^Sb. In the spectrum from 
Rev. 0356 taken during flare, the Sun and Earth atmosphere 
emission by light isotopes is observed at high energy. There 
are five lines at about 2200, 3400, 3900, 4150, and 4450 keV, 
which can be identified with the emission from ('^O, "B, "^B, 
"■O, i^N, "N), ('^N), ('^C, ' ^N), ("'O) and ('^C, ^B), re- 
spectively (^S hare & MurphvlbOOl ) . During quiet Sun states, the 
Earth atmospheric emission is not observed, as demonstrated by 
analysing the spectra extracted from Revs. 0401, 0404, 0404, 
and 0406, when the IN TEGRAL Earth occultation observation 
(IChurazov et al.ll2007l) was performed. 

3.2. Energy calibration 

A routine energy calibration for PICsIT is based on the special 
data (so-called S5 datastream) collected in coinciden ce with the 
emission from the onboard ^-Na calibration source (Bi rd et all 
2003). S5 data are collected as sep arate 64-channel spectr a and 
monitored by the instrument team jMalaguti et alj|2003bh . pro- 
viding the parameters for the channel-energy conversion. This 
standard calibration assumes that energy is a linear function of 
the channel number, which in general is not necessarily true. 
Studies of 256-channel background spectra completed directly 
by using the standard spectral-imaging histograms of single 
events can help in controlling the energy calibration, as well as 
extending it outside the 51 1-1274 keV range of S5 data. 

Energy calibration tests were completed as follows. 
Background spectra were extracted for all public data from revo- 



p. Lubinski: Analysis of INTEGRAL/PIC sIT data 



5 



1 1 1 1 1 1 


1 1 1 1 1 1 
0052 (normal) 






, — , — \ — 


-H > > > \— 

0130 (after flare) 






, — , — 1 1 1 1 

■ \ 


-H < ' ' h- 

0356 (flare) 


, — , — \ — 

■ + 

_ t 




-H ' ' ' \— 

0541 (normal) 

.1 . . . 1 



500 1000 5000 

Energy [keV] 



Fig. 3. PICsIT background spectra from Revs. 0052, 0130, 0356, 
and 0541. Data are divided by the power-law model with a spec- 
tral index T set to 1 and a normalization at 1 keV equal to 1 .25 
keV ' cm"^ s Spectra representing non-flaring Sun periods 
(0052, 0541) are separated by exactly 4 years. The spectrum 
from Rev. 0130 was taken 3-5 days after the most powe rful so- 
lar fla re ever recorded (X28-I-, 4 Nov 2003, Sp aceWeather.comI 
I2007h . when the Sun already returned to quiescence. The Rev. 
0356 spectrum shows the background emission observed 5- 
7 days after the fifth l argest solar flai-e (X17, 7 Sep 2005, 
ISpaceWeather.coml2007l) . while there was still considerable flar- 
ing activity. 



lutions 0039-0676. The positions of the 511, 667, and 1 779 keV 
lines were then fitted in channel space using Xspec 1 1 (lArnaudI 
11996 ). The lines were modelled with Gaussian shape and the 
local continuum was asummed to be in the form of a power- 
law. After some initial calibration tests, the most reliable gain 
value was established to equal the width of one basic channel, 
i.e., 7.1 keV, in agreement w ith the standard ener gy calibration 
performed by the IBIS team dFoschini et al.ll2007l) . The only re- 
maining parameter needing to be adjusted was then an offset, a 
shift in energy units to be applied to the original onboard his- 
tograms to ensure that the three tested lines were in the correct 
positions. The best-fit offset was found to be -1-9.6 keV with re- 
spect to the standard energy calibration of PICsIT. The first us- 
able channel of spectral-imaging single-events histograms was 
channel number 10, which in the standard OSA 7.0 calibration 
starts at 203 keV and in the new calibration presented here at 
212.6 keV. 

The mean positions and 1-cr errors found for the background 
lines in the spectra from revolutions 0170-0676 are equal to 



509.2+5.0 keV, 668.5+3.8 keV and 1772.8+4.1 keV. Before rev- 
olution 170, a different binning was applied to the onboard his- 
tograms and the data were treated separately. The corresponding 
mean positions of the calibration lines were found to be shifted 
by about +5 keV with respect to the results quoted above. 

The main factor driving the offset of energy calibration i s 
the temperature of the PICsIT detector (Malag uti et al.l l2003b). 
Short time variations in the line position reach as high as about 
15 keV on timescales of weeks. This is a large value compared 
to the minimal width of the bin used in spectral extraction, 21.3 
keV. In principle, this can introduce fluctuations in the extracted 
count rates. As we show in Sect. 16.31 Crab count rates are stable 
over the mission time, which indicates that the offset correction 
is not very important, although it should be taken into account in 
a proper preparation of the background maps. This is currently 
achieved by preparing several maps for revolutions with large 
temperature variations, using data from periods with an approx- 
imately constant temperature. 

3.3. Variability 

The count rate of the PICsIT background varies strongly in re- 
sponse to the cosmic-ray intensity changes, since it is moder- 
ated by the decay time of various radioactive nuclei produced 
in the detector neighborhood. On average, the background count 
rate increased with time during the first four and a half years 
of INTEGRAL operation, as shown in Fig.lD This increase was 
caused by both the solar modulation of the cosmic -ray intensity 
and the accumulation of the radioactivity in the detector assem- 
bly. The increase in the background count rate is irregular, there 
are periods of high variations but there are also periods of almost 
constant count rate. Since March 2007 (around Rev. 0540), the 
background appeared to have saturated at a level of about 60% 
higher than the initial one. However, below 250 keV the increase 
in the background count rate is far more significant, higher than 
100% since the beginning of the mission and appears still not to 
be saturated. 

Within a single revolution, a substantial variability is also 
observed, induced mainly by the spacecraft passages through the 
Earth radiation belts and the Solar activity. An example of this 
variability is presented in Fig. |5] for six energy bands between 
213 keV and 1940 keV. This figure illustrates a rather typical 
behaviour, where the total count rate decreases with time for en- 
ergies below w511 keV and increases above that energy. The 
very low energy bands (close to 200 keV) usually exhibit a dif- 
ferent type of variation, due to the track events. The second and 
third (from top) panels of Fig. |5] indicate that the background 
variability can be quite different in two close energy bands. 
Consequently, spectral extraction with a background model con- 
structed for too wide an energy band can be inefficient when the 
radiation in the low-energy part of the band increases but in the 
high-energy part decreases with time. In general, there is a vari- 
ety of background light curves observed, sometimes the changes 
in all energy bands are correlated, sometimes there are flare-like 
events or rapid count-rate variations that are only observable for 
a limited energy range. Therefore, for a correct background mod- 
elling, the variability in narrow energy bands must be determined 
for each revolution. 

3.4. Background maps 

The background model for PICsIT is the map of pixel count 
rates corresponding to the background emission registered in a 
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Fig. 4. Variation in the PICsIT background count rate during the 
mission. Error bars correspond to the standard deviation of the 
mean for a given revolution and are much larger during solar 
flaring periods. 



given energy band during some period of time. Both the sum 
of photons originating in all point sources in the field of view 
and the summed difiiise sky emission equal at most of the or- 
der of several per cent of the total detector count rate at energies 
above 200 keV. Accordingly, the background map consists of a 
sum of shadowgrams collected during a number of dithering ob- 
servations, such that the point sources and the diffuse emission 
patterns are spread over the detector surface. In this way, the 
background model contains the same data used subsequently to 
extract the source count rate. The dithering observational strat- 
egy thus allows us to obtain for the coded-mask instruments 
the highest quality background estimate, based just on the ac- 
tual data. This is one of the most important virtues of masked 
pixellated detectors compared with standard telescopes, where 
the background must be determined in some indirect way, such 
as a rocking observing strategy as applied for OSSE. 

While there is no variation in count rate within each pixel 
during an observation, the background map provides the most 
accurate, experimentally-determined background model, with an 
accuracy increasing with the exposure time. This situation never 
occurs in reality and the main task in achieving a correct back- 
ground model is to describe properly its variation in time. Some 
observations are also taken in special conditions, such as star- 
ing mode, when the map cannot be prepared with the contempo- 
rary data. Consequently, the variability model should allow us to 
check if another map can be used for these special observations. 
The simplest way of taking into account the background vari- 
ability in time is to assume that there is only one global trend, 
coherent for all pixels. The only parameter that has then to be 
determined with the source count rate is the background map 
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Fig. 5. Total PICsIT count rate measured in eight energy bands 
during Rev. 0159. IJD is the INTEGRAL Julian Date (IJD = 
MJD -51544.0). 



normaUzation for a given science window. As shown in Sect. |6] 
this approach produces completely acceptable results. 

In view of the considerations mentioned above, the only pa- 
rameter that is free in the modelling of the background is the 
time period for which the background map is prepared. Since 
the spacecraft activation is caused predominantly by passages 
through the Earth radiation belts, the natural choice of time inter- 
val for which a map is prepared is the revolution period. Owing 
to the large number of counts collected by each pixel during this 
time period, the accuracy of the revolution-averaged background 
rate is sufficiently high. On the other hand, if there is a sub-orbit 
variation observed, the shortest time for which the background 
map can be prepared is limited by the condition of smoothing 
the source illumination pattern. This will depend on the actual 
dithering pattern, but usually about 10 science windows are suf- 
ficient to obtain a sufficiently smoothed summed shadowgram. 
Nonetheless, producing more than one background map for a 
single revolution is unnecessary, unless there is clear variability 
in the background mean count rate or the detector temperature 
measured during the orbit time. In practice, almost all results 
presented later were obtained with the use of a single map, apart 
from several observations of Crab and Cyg X-1, where two maps 
per revolution were needed to suppress an artificial trend found 
in the source light curves. 

Data used to prepare the map were selected according to the 
mean count-rate level over all pixels in a given science window. 
The fluctuations in that value usually did not exceed 3% for pe- 
riods with a stable background and this criterion was adopted 
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for the PICsIT data analysis presented in this paper. In the spec- 
tral extraction, a similar criterion was used, excluding the data, 
when the background map normalization was outside the 0.97- 
1.03 range. This selection automatically excludes periods with 
solar flares or science windows affected by the passage through 
the radiation belts at the beginning or end of revolution. 




Fig. 6. Distributions of differences between 336-448 keV shad- 
owgram from Rev. 0079 and corresponding shadowgrams from 
Revs. 0070 to 0078 and 0080 to 0090 (starting from the upper 
left, row by row). Warmer colours (yellow, red) correspond to 
the higher count rates. The smallest differences (flattest images) 
are observed for Revs. 0072, 0075, 0078, 0080 and 0081. There 
was a small solar flare in Rev. 0076 that changed also the back- 
ground distribution in the following revolution. 

When it was not possible to prepare the map for part or an 
entire revolution, one could use, in principle, a map prepared for 
another period close in time. Tests done for the Crab indicated 
that this approach did not always provide good results, especially 
when the map is constructed with data from another revolution. 
The simplest way of confirming the adequacy of a map is an in- 
vestigation of the differences or ratios between the pixel count 
rates in the summed shadowgrams from both the map and the 
source observation. Figure |6] presents an example of this test 
made for the summed shadowgram from Rev. 0079 when Cyg 
X-1 was observed. The flattest difference shadowgrams are those 
for Revs. 0078 and 0080 (two rightmost in the second row of the 
plot). It has to be kept in mind, however, that this test may not 
be applicable when a relatively strong source is observed in star- 
ing mode, because then there will be a systematic imprint on the 
shadowgram from the source illuminating the same pixels all of 
the time. A possible solution to this limitation is to test the back- 
ground stability over several revolutions adjacent to the actual 
source observation and to prepare some mean map. 

4. Count-rate extraction method 

The necessity of merging a large amount of PICsIT data to detect 
weaker sources not only compels one to model the background 
precisely but also to apply a method that is strictly correct for 
data with a low signal-to-noise ratio. Otherwise, the inadequacy 
of the method can produce results with a bias comparable to the 



signal strength. The most common way of estimating some pa- 
rameter value and its confidence limits is the application of the 
test statistic. Although it is well known that the conditions 
needed to justify the use of the test are not fulfilled for low 
number of counts, it is worth reminding us of all limitations of 
the standard technique. 

1. The assumption that the Poisson distribution can be ap- 
proximated by a Gaussian distribution for a large number of 
counts is not always justified because these distributions differ: 
the Poisson distribution is asymmetric with non-zero skewness 
and kurtosis, whereas the Gaussian distribution is symmetric 
with zero skewness and kurtosis. When the signal has an ampli- 
tude comparable with the difference in shapes of the Poisson and 
Gaussian distributions for the total (signal +background) number 
of counts, the Gaussian approximation will provide a systemati- 
cally biased result. 

2. The Poisson distribution is defined for parameters of non- 
negative values. Count rates also cannot be negative and any 
analysis using a Gaussian or other distribution in negative sig- 
nal space does not offer a completely proper description of the 
phenomenon. 

3. The difference of two parameters both characterized by a 
Poisson distribution is not Poisson distributed. Consequently, the 
signal+background and the background alone should be mod- 
elled simultaneously instead of subtracting the background rate 
from the total rate, especially when this subtraction leads to neg- 
ative signals. 

4. Since the true variance is usually unknown, the observed 
number of counts serves as a variance estimate for the statis- 
tic. This approximation is invalid for a small number of counts 
and many solutions have been proposed to improve the variance 
estimate, but none of them performs well for weak signals, when 
zero net counts are measured in a majority of trials. 

5. To estimate the uncertainty of a parameter of interest when 
many parameters are fitted simultaneously, one usually consid- 
ers the limits corresponding to th e change in the statistic by a 
certain value (iLampton et al.|[T976i) . This approach based on the 
projection of the respective parameter confidence region onto the 
parameter axis often overestimates the error. Moreover, its basic 
condition, namely the use of the model variance, is not met in 
many practical applications, as stated in point 4. 

The solution to several issues listed above is application of 
a statistic proposed by Cash ( 1979) based on a maximum like- 
lihood estimation for the Poisson distribution. The application 
of the Cash statistic allows us to avoid problems arising due 
to the differences between Poisson and normal distributions but 
still has some limitations that will be discussed in Appendix A. 
Due to its direct handling of probability density function (here- 
after PDF) associated with the model parameters, a Bayesian ap- 
proach provides a more general and efficient solution to many 
statistical questions, including the 'low-number-of-counts' is- 
sue. 

A Bayesian technique for dealing with the Poissonian type 
data in th e presence of k nown or unknown background was pre- 
sented b v iLoredol (11990). This type of approach applied to a low- 
quality ASCA data obtained for weak AGNs, was demonstrated 
to re construct preci sely the assumed shape of the iron line com- 
plex (lLubinskill200 4). Due to its efficiency, the method has now 
been generalized to deal with the extraction of count rates for ob- 
jects observed with PICsIT. Developments have occured in two 
respects: the Poisson PDF is now computed in two (or more) 
dimensions, and the model now takes into account the main ob- 
servation conditions. 
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Table 2. Typical source (R,) and background (/?;,) count rates 
observed with PICsIT (upper part. Rev. 0301) and ISGRI (lower 
part. Rev. 0079) for sources of 50 mCrab and 2 mCrab strength, 
respectively. Corresponding numbers of counts, A^, and Ni,, were 
calculated for science windows with a 2 ks exposure time. S/N is 
the signal-to-noise ratio, and kio- is the number of science win- 
dows with the 2 ks exposure time needed to obtain a 3-sigma 
detection. 



Energy range R, N, Rt Ni, S/N k^^ 
[keV] [cps] [cps] 



The probability density function puis, bk) of the source count 
rate parameter s and the background map normalization factor 
bk in the science window number k for a given energy band is 
defined as the product of Poisson PDFs P for all tik active pixels 

Pk{s, bk)^—Y\ PiNdi); (m^di) + bkB(i))Mi)), (1) 

where the factor Q normalizes the pk{s,bk) integral to 1, Nkif) 
is the number of counts measured by the /-th pixel, rjkii) is the 
source PIF value for that pixel, corrects s for the off-axis ef- 
fect, and B{i) is the efficiency-corrected background count rate 
given by the background map. The factor fkif) converts count 
rates into counts, 

fk{i) = Tks,{i)l{32 X 64), (2) 

using the exposure time Tk and the pixel efficiency et(/), where 
the source count rate is normalized to the half of the detector 
area (32x64 pixels). 

The density pk{s, bk) is the posterior probability density, i.e. 
the prior probability densities on s and bk are assumed to be 
uniform. When a large amount of data is used, as in the PICsIT 
case, with thousands of pixels and hundreds of science windows, 
the final density distribution associated with the detected source 
signal is sufficiently concentrated not to be affected by the choice 
of the prior distribution. 

Since the calculation of the PDF must be performed for a 
relatively wide range of s and bk, the product in Eq. [T]is com- 
puted by summing the PDF logarithm to avoid the multiplication 
of many small numbers. The final PDF p{s) is derived after in- 
tegrating pk{s, bk) over the background normalization parameter 
and then computing the product of pk{s) for all science windows 
of interest. Using p(s), one can extract all needed parameters 
such as the mean, the median, or the desired credibility inter- 
vals. For PICsIT, it is possible to have a detection in a single 
science window only for sources at least as bright as the Crab. 
To achieve a reliable result for weaker objects, a large number 
of science windows must be merged. Table |2] compares typical 
values of the number of counts, count rates, and signal-to-noise 
ratios for data of a 50 mCrab source observed with PICsIT and 
a 2 mCrab source observed with ISGRI. 

Examples of Poisson PDFs obtained for a true source are pre- 
sented in Fig.|7] PDFs computed for a single science window are 



very broad when compared to the source count rate. For an ob- 
ject as weak as GRS 1758-258 (about 70 mCrab in the 300-500 
keV band), even one-revolution of data is not enough to achieve 
a detection. Nevertheless, thanks to a frequent monitoring of the 
Galactic centre field by INTEGRAL, it is possible to follow the 
seasonal (i.e. on the scale of several months) changes in the high 
energy emission from this source in several energy bands below 
500 keV. 
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Fig. 7. Examples of Poisson PDFs obtained for GRS 1758-258 
in the 298-319 keV band. Left: broad distributions obtained for 
several single science windows from Rev. 0105 (dashed lines) 
and a narrower PDF from all Rev. 0105 data merged (solid line). 
Right: The same total PDF from Rev. 0105 (dashed line) com- 
pared to that of Rev. 0103 (dot-dashed line), and the final result 
for all data of August - October 2003 period (solid line), with an 
exposure time of 1.95 Ms. 

The technique presented in this paper is the most cor- 
rect means of handling the Poisson-distributed data within the 
Bayesian framework. Due to this first application of such a 
method to extract the astrophysical y-ray data, some tests of the 
performance are needed to check for example the limiting signal 
that can be extracted and compare this technique with several 
other commonly used methods. An obvious test is to use this 
method for the ISGRI spectral extraction and to compare the re- 
sults with those of the standard OSA software. The results of 
these tests are presented in Appendix A. 

The computation time needed to extract the source count 
rate over a PDF grid as fine as possible (the grid step should 
be much smaller than the expected rate) is long compared to the 
time needed by the method using the fitting with some statis- 
tic. For only two parameters handled (the count rate for a sin- 
gle source and the background normalization), computing over 
a PDF grid can be performed even for a large data set. However, 
if there are several sources in the field of view, the computation 
over the grid has to be replaced by some Monte Carlo integration 
of multi-dimensional PDFs. Such procedures are time consum- 
ing and need some adjustments to accelerate the data analysis. 
Due to the small number of sources observed with PICsIT, this 
type of software does not appear to be needed, especially in a 
situation when the emission from any contaminating source is 
far weaker than the background. Nevertheless, a version of the 
software with the Monte Carlo integration implemented was pre- 
pared and tested with the simulated data for up to 7 sources in the 
field of view. These tests show that only in the case of sources 
separated by a short distance (less than several degrees) there can 
be some contamination for weaker sources. An independent test, 
where the single-source PDF method was applied separately to 
ISGRI data of NGC 4151 and NGC 4051 (separated by 5.3°), 
has shown that the results for both objects are fully consistent 
with those of the standard OSA software handling many sources 
at once. 
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5. Detection limits 

The standard way of verifying whether the observed signal cor- 
responds to a detection is the 'n-sigma' significance test. The 
test statistics (significance) is defined as the difference between 
the measured signal+background and the background estimated 
in an independent way, divided by the noise represented by the 
standard deviation of that difference. Detection can be claimed 
if the observed significance is higher than the threshold for a 
given probability level, e.g., 3 for probability level 0.997 (3- 
sigma test). The significance level is often set to 5- or 6-sigma 
because such a conservative approach is supposed to account for 
the number of trials or to balance all possible systematic effects 
not incorporated into the uncertainty calculation. 

In true situations, the significance estimate should take into 
account all circumstances of the measurement: the detector 
model, the background model, and the method used to extract 
the count rate. For coded-mask instruments, one has to take into 
account the mask pattern and transparency, the spatial resolution 
of the detector and the background non-uniformity. There are 
many approaches developed to determi ne analytically the sensi- 
tivity of coded-mask instruments (see ISkinneJbOOSl and refer- 
ences therein). 

This simple 'n-sigma' significance test cannot be applied to 
the results of the technique presented here. The method based 
on the Poisson PDFs extracts the net source count rates with- 
out direct background subtraction and operates in a physical, 
non-negative parameter space. The mean net count rate deter- 
mined by integrating a PDF over a non-negative argument can- 
not be zero even when there is no signal, although it decreases 
with increasing exposure time. Therefore, there will always be 
some signal excess above zero resulting in an overestimation of 
the true significance level. On the other hand, for the weak sig- 
nals expected to be observed with PICsIT, a too conservative 
approach might exclude some potential sources from the analy- 
sis. Therefore, we would like to find a way of estimating the true 
noise level. This will be useful not only for the verification of the 
results but will also help to judge if the source with a given flux 
level can be detected after a given observation time. 

The question addressed here is: "what is the probability of 
obtaining a false signal, arising from the background (noise) 
fluctuation and described by the probability density function 
with a given mean count rate, for an observation lasting a given 
time". To find the answer, we compiled the distribution of mean 
count rates extracted for random positions in fields that do not 
include strong sources. The tail of the distribution was then fit- 
ted to find the upper noise limit corresponding to a given prob- 
ability level. Tests of the PICsIT noise level were completed for 
six wide energy bands: 277-362, 362-461, 461-632, 632-930, 
930-1938, and 1938-3131 keV. The limits of these wide bands 
were adjusted to have a coherent background variation within 
a given band (see Sect. 13. 3t . Input empty-field data were se- 
lected as observations when all of the eight possibly brightest 
PICsIT sources (Crab, Cyg X-1, GRS 1758-258, 4U 1700-377, 
GRS 19I5-H105, Cen A, IE 1740.7-294, XTE J1550-564) were 
farther than 45° away from the pointing direction. This criterion 
provided 5599 science windows from Revs. 0047-0450 with the 
total exposure time equal to 14.48 Ms. Input PIF data were cal- 
culated for 91 (13x7) possible positions corresponding to the 
hexagonal dit hering pattern w ith a so-called "wandering cen- 
tre of pattern" ('K uulkersl2005h . The hexagonal dithering pattern 
was chosen because in this case the object is in the fully coded 
FOV all the time and there is no need to correct the exposure 
time for incomplete detector illumination. 



At the beginning of each test, the PIF was randomly selected 
for each input science window. The net source count rate was 
then extracted for each science window using the method de- 
scribed in the previous section. These results were merged for a 
subset of all science windows to obtain the mean count rate for a 
given assumed total exposure time. For shorter times, it was pos- 
sible to select the subset several times because the entire sample 
was sufficiently large to avoid repetition. To ensure that the tail 
of the mean count-rate distributions could be determined to a 
good precision, approximately 1000 tests had to be performed. 
However, since computation of one test takes several hours on 
a 20-CPU grid network, it is impossible in practice to complete 
1000 tests for times longer than 1 Ms. 

Upper noise limits were calculated for probability values 
corresponding to the commonly used 'n-sigma' significance lev- 
els, where n was in the range 1-5. The limit of 1-sigma signif- 
icance level was determined as the mean count rate value for 
which the sum of runs, starting from the minimal mean count 
rate, reaches 68.3 % of the total run number. A model function 
was then fitted to the tail of the distribution, i.e. above a 1-sigma 
limit and the 2-5 sigma limits are found via integration of that 
tail model. Because the true shape of the distribution was un- 
known, three model distributions were tested for the 277-362 
keV band, where the number of tests was much larger, reaching 
more than 1 1 000 runs for the shortest 10 ks exposure time. The 
Gaussian function was used to model a rapid decline and the ex- 
ponential function to model a slow decrease. The third model in- 
cluded a gamma distribution, which reproduced an intermediate 
decrease. For distributions analyzed with a large number of tests 
(> 1000), the best fit was always obtained with the Gaussian dis- 
tribution. The fitted gamma distribution was usually quite similar 
to that of the Gaussian, but the 3-sigma upper limit was found to 
be slightly (1-5%) larger The exponential model, in turn, seems 
to underestimate the decrease in the tail, giving always the high- 
est The uncertainty in the fitted 3-sigma upper limit was es- 
timated to be, in general, not larger than 10% and 20 % for the 
exposure times below and above 1 Ms, respectively. This esti- 
mate was based on tests with an increasing number of runs using 
different tail models. 

Figure[8]presents the example of noise limit estimate for the 
277-362 keV band and exposure times of 10 ks and 1 Ms. The 
distributions of mean count rates extracted for empty fields are 
asymmetric and clearly peaked above zero. Comparison of the 
three fitted tail models is shown in the upper panel of the figure. 
For 1 Ms, three PDFs are presented with the highest mean count 
rate, which illustrate the typical shape of the probability density 
distribution at this count rate level. These PDFs are compared 
with two PDFs obtained for the microquasar GRS 1758-258 af- 
ter observations lasting about 1.05 Ms. The GRS 1785-258 PDFs 
correspond roughly to 5 and 6-sigma detections. 

The results of the noise level studies for 3-sigma detection 
limits are collected in Table [3] A wide range of effective expo- 
sure times should allow for a relatively precise interpolation for 
a particular observation time. The last row of the table shows the 
mean PICsIT fluxes measured for the Crab that can be used to 
convert the limits into Crab units. PICsIT sensitivity limits for 1 
and 10 Ms are compared in Fig. |9] with the 1 Ms sensitivity lim- 
its computed in an analytical way for PICsIT, SPI, and OSSE. 
Sensitivity limits for PICsIT and SPI were calculated using the 
on-line Observation Time Est imator (OTEfl The formulae used 
by the OTE were described in Belanger (2008). OSSE sensitivity 
data are taken from Fig. 3 of .Winkler, (, 1 994) . 
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Fig. 8. Distributions of mean count rates in the 277-362 keV 
band extracted for fake sources in empty fields. Upper panel: 
Total exposure time 10 ks, 11 000 tests. Solid, long-dashed and 
short-dashed lines show the Gaussian, gamma, and exponential 
models, respectively, fitted to the tail of the distribution, i.e. 
above the count rate corresponding to a 1-sigma upper limit. 
Vertical lines indicate the positions of the best-fit 2 and 3-sigma 
upper limits. Lower panel: Total exposure time 1 Ms, 1000 tests. 
The solid line shows the Gaussian function fitted to the tail of the 
distribution. Vertical lines present the 1 to 5 sigma upper noise 
limits computed using the fitted Gaussian tail. Three extremely 
high-noise distributions for 1 Ms are shown with dotted lines, 
and two examples of count rate distributions obtained for GRS 
1758-258 with a similar observing time are drawn with dashed 
lines. 

Table 3. PlCslT 3-sigma sensitivity limits determined through 
the extraction of count rates for fake sources in the empty field 
observations. The last row presents mean fluxes determined for 
Crab. 



Eff. exp. Energy band [keV] 
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Despite the different approach used in the computation, 
PICslT sensitivities based on the noise level tests and given by 
the OTE tool are consistent for the first two energy bands. Above 
460 keV, the OTE estimates are well below those determined ex- 
perimentally. Since noise determination based on the empty field 



observations takes into account several effects not included in the 
OTE model (e.g. background non-uniformity, background level 
evolution with time, and a realistic PIF model) one can expect 
OTE to overestimate the PlCsIT sensitivity. The same concern 
can possibly be raised for the SPI sensitivity estimates made in 
an analytical way. The PlCslT detector should be more sensitive 
to continuum observations because of its larger area and vol- 
ume, the higher quantum efficiency of Csl compared with Ge 
crystals in terms of registering photons, and higher number of 
detector and mask pixels allowing for a more robust disentan- 
glement between the signal and the noise. Nevertheless, a more 
reliable comparison of the SPI and PICsIT sensitivity is not pos- 
sible unless a noise study similar to that completed for PlCslT 
is also completed for SPI. On the other hand, this comparison is 
probably more secure for the OSSE detector because of the rock- 
ing strategy of observations with that instrument. PlCslT appears 
to be much less sensitive than OSSE. As can be seen in Fig.|9] 
PlCslT needs about 10 times longer exposure time to reach the 
OSSE S/N ratio because its volume is three times smaller This 
lower efficiency is partly compensated, however, by a larger frac- 
tion of observing time (because of a higher orbit). 
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Energy [MeV] 

Fig. 9. PICsIT sensitivity limits for a 3-sigma detection for an 
exposure time of 1 Ms. The thick solid lines show the estimates 
based on the noise level test for 1 and 10 Ms effective exposure 
times. Dotted and dash-dotted lines correspond to the OTE tool 
results for 1 Ms observation with SPI and PICsIT, respectively. 
The dashed line shows the 1 Ms sensitivity limits for OSSE. 
Thin lines show the fluxes for 10 mCrab, 100 mCrab and 1 Crab 
sources with the spectral slope fixed to 2.2. 



6. PICsIT spectra 

We now review the PICsIT spectra extracted for various objects. 
Some of these spectra are used to verify PICsIT calibration, by 
comparing them with the spectra obtained from two other high- 
energy INTEGRAL detectors, SPI, and ISGRI. Potential prob- 
lems that might arise due to systematic effects are checked us- 
ing Crab spectra. The results of the spectral extraction method 
based on the Poisson PDF (hereafter PPDF) technique described 
in Sect. |4] are compared with the results obtained with the stan- 
dard OSA 7.0 software. Finally, examples of spectra derived for 
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several detected objects are shown against the noise level ex- 
pected for a given exposure time. 

6.1. Calibration 

Cross-calibration tests are one of the most important elements 
of the verification of a given detector performance. Due to the 
presence of three high-energy detectors onboard INTEGRAL 
observing the source at the same time, these tests are easier 
and more reliable than in the case of satellites hosting only a 
single instrument operating in a given energy range. Moreover, 
the SPI detec tor was carefully calibrated before the launch of 
INTEGRAL ( iAttie et alj|2003l) . which ensured that the calibra- 
tion check was independent of the details of the particular model 
assumed for a calibration source such as the Crab. A lot of effort 
has been invested since the beginning of the mission to calibrate 
well the ISGRJ and SPI instruments, resul ting in a good perfor - 
mance in the energy range below 100 keV dJourdain et al.l2068l) . 
Due to a limited ISGRI sensitivity above 100 keV, there is how- 
ever still a spectral slope issue remaining for the high energy part 
of the spectra. On the other hand, PICsIT standard OSA spec- 
tral extraction was applied only to a limited number of sources. 
This enabled only a crude check of the calibration, suggesting 
that there was a n overall agreement b etween PICsIT and ISGRI 
spectral results ('F oschini et al.ll2007l) . 

The first PICsIT response files for single-event spec- 
tra were provided to the user community with the OSA 
4.2 release in December 2004. The two response files 
currently used are the RMF pics_srmf_grp_0003.fits (or 
pics_srmf_grpj)005.fits rebinned to the default energy bins) and 
the ARE pics^arfj:sp_0003.fits. Both files were used in the tests 
presented here, although both interpolation and smoothing were 
applied to them to correct for discontinuities caused by the 
too wide energy bands used in their Monte Carlo computation. 
These corrections were necessary for obtaining smooth back- 
ground spectra without jumps, obse rvable when the standard re- 
sponses are used (see e.g. Fig. 6 in lBelange3l2008l) . The source 
spectra with finer binning are also smoother after this correction, 
although the overall calibration remains similar to the standard 
one. The final change that remains to be applied to the current 
response is the correction for the influence of the mask pattern. 
Preliminary tests showed that this correction affects mainly the 
high-energy part of the response, with the spectral fitting results 
below 1 MeV remaining virtually the same (L. Natalucci, private 
communication). 

Another issue related to the INTEGRAL cross-calibration 
is the ISGRI spectral slope above 100 keV. Due to superior 
sensitivity up to about 150 keV, ISGRI spectra dominate the 
INTEGRAL spectral fitting. Therefore, even a small discrep- 
ancy between the spectral slope of data from ISGRJ and the two 
other detectors, can change substantially the results of a broad- 
band fit, producing a model with an incorrect spectral shape at 
high energy and incorrect relative normalization between instru- 
ments. For a long time, the photon index F of the reference Crab 
model used in ISGRJ calibration was set to be 2.225 for the 
entire ISGRI energy range. In OSA 7.0, following the recom- 
mendation of the INTEGRAL Users Group, it was changed to 
2.1 below 100 keV and 2.34 above 100 keV, to match the SPI 
Crab spectrum, as provided by the SPI team. However, the latter 
value is clearly higher than F « 2.2 inferred from the Crab spec- 
tral analysis with the standard OSA SPI software, provided in its 
version 7.0 and several earlier releases. The l atest INTEGRAL 
cross-calibration report (iJourdain et al.l l2008h presented a new 
result for SPI: with the updated spectral extraction software used 



Table 4. Cross-calibration between the INTEGRAL high-energy 
instruments. A power-law model was fitted to the Crab (three 
upper rows) and Cen A (three lower rows) spectra. The best- 
fit parameters for Crab are: a photon index F of 2.19+0.01 and 
a normalization at 1 keV of 17.1+0.5 photons keV"' cm"^ s"' 
with a reduced of 1.71 for 49 degrees of freedom. The best- 
fit parameters for Cen A: a photon index F of 1.80+0.01 and 
a normalization at 1 keV of 0.153+0.006 photons keV ' cm"^ 
s ' with a reduced of 1.16 for 61 degrees of freedom. No 
systematic errors were added to the statistical ones. 



Instrument 


Exposure 


Energy range 


Rel. norm. 




[ks] 


[keV] 




ISGRI 


1168 


100-462 


0.875+0.003 


SPI 


910 


105-957 


(1.00) 


PICsIT 


981 


277-930 


0.75+0.01 


ISGRI 


771 


19^62 


0.89+0.02 


SPI 


584 


26-352 


(1.00) 


PICsIT 


429 


277-454 


0.85+0.13 



by the SPI team, the Crab spectral slope above 100 keV is now 
equal to 2.22, i.e. a value fully consistent with the OSA result. 
One can thus expect the next OSA releases to provide revised 
ISGRJ response files, adjusted accordingly. To anticipate this 
inevitable result, the standard ISGRI OSA 7.0 ARFs were cor- 
rected to match the value 2.22 above 100 keV for the purpose of 
the testing described below. 

Calibration tests were made for two sources. Crab and Cen 
A, the latter to enable checks to be completed for a weaker ob- 
ject and with a different spectral shape. High energy spectra were 
produced with a large data set to achieve the highest possible 
precision. For Crab, all data from the standard dithering observa- 
tions made in Revs. 0239, 0300, 0365, 0422, 0483, 0541, 0605, 
0665, and 0666 were used. This choice takes into account the 
fact that before Rev. 0216 SPI had not yet lost two of its detectors 
and its response was thus different during earlier observations. 
In the case of Cen A, it was impossible to obtain acceptable SPI 
spectra for observations completed between January 2005 and 
February 2007. Unfortunately, in this period there were strong 
background variations caused by enhanced Sun activity. PICsIT 
spectra for that period were also clearly worse than those ob- 
tained for earlier observations. The Cen A calibration spectra 
comprise therefore only data from observations made before 15 
February 2004 (up to Rev. 0163), for which the source was de- 
tected at energies up to about 460 keV. Calibration tests with 
ISGRJ and SPI spectra of the Crab were limited to the energy 
band above 100 keV, to decrease the influence of the low energy 
part on the fit completed for the PICsIT spectral band. 

Results of spectral fitting performed with the calibration 
spectra are presented in Table |4] and illustrated in Figs. 1101 and 
nn In agreement with the results of tests done with the stan- 
dard OSA software, the PICSIT calibration appears to be re- 
liable, validating the instrument model included in the Monte 
Carlo software used for the response simulations. Both the Crab 
and Cen A spectra measured by three INTEGRAL instruments 
are fully consistent, producing similar spectral slopes also when 
fitted separately. The cross-calibration factors measured with re- 
spect to SPI obtained for ISGRI and PICsIT are both below 1. 
The PICsIT cross-calibration factor is poorly constrained for the 
Cen A spectra, where the SPI and PICsIT useful bandwidths 
are rather narrow. The relative PICsIT/ISGRI normalization ob- 
tained for the Crab of 0.86 is well above the value of 0.52 quoted 
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for PICsIT spectra e xtracted with the standard OSA software 
dFoschini et al.l2007l) . This issue is discussed further in Sect. 16. 31 



1 r 
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Fig. 10. Cross-calibration test completed for the INTEGRAL 
Crab spectra. ISGRI and SPI spectra were extracted with the 
standard OSA 7.0 software, and the PlCSlT spectrum was ex- 
tracted with the PPDF method. Results of a power-law fit are 
given in Table |4] 
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Fig. 11. Cross-calibration test for the INTEGRAL Cen A spec- 
tra. ISGRI and SPI spectra were extracted with the standard OSA 
7.0 software, and the PICSIT spectrum was extracted with the 
PPDF method. Results of a power-law fit are given in Table H] 



6.2. Other tests 

Besides the verification of the PICsIT calibration, other tests 
were needed to check the performance of the spectral extraction 
software related to the instrument model. One of the main effects 
to be tested was the dependence of the results on the off-axis an- 



gle at which the source was observed. Many possible problems 
can arise from different orientations of the spacecraft with re- 
spect to the object. It is impossible to account for all of them in 
the instrument model included in the spectral extraction or in the 
response-simulation software. Problems with the off-axis obser- 
vation results were reported for the standard OSA software and 
some corrections to the sky-image deconvolution code were in- 
troduced starting from the OSA 6.0 release (Foschini 20071^. 
The PIF-based spectral extraction presented in this paper ac- 
counts for some off-axis effects, as akeady mentioned in Sect. 
|2] The Crab spectra presented in Fig.[T2]show that these correc- 
tions are sufficient to ensure stability of the results. They were 
extracted for all data taken during and after Rev. 0170, and for 
several subsets of those data, selected according to the off-axis 
angle. Both large off-axis angle (above 10°) and small off-axis 
angle (in the 0.5°-l° range) spectra are fully consistent with the 
total spectrum. However, it appears that the spectra produced 
from the on-axis data (off-axis angle below x 0.3°) are quite dif- 
ferent. They are found to exhibit strong spectral hardening ob- 
served in this case, starting at an energy of w 400 keV and re- 
sulting in a Crab spectral index F of around 1 .4. A similar effect 
is observed for on-axis spectra extracted with the standard OSA 
7.0 software. 

Deeper studies of the on-axis phenomenon indicate that its 
amplitude in single-revolution spectra is variable. The harden- 
ing is sometimes strong or very strong, but sometimes it is quite 
weak or not present at all. In addition, the hardening appears at 
different energies for individual spectra. A possible explanation 
of at least part of this effect is discussed below, but a full ex- 
planation will need more tests and additional data. For instance, 
a micro-dithering observation could be useful. Nevertheless, in 
practice, sources other than the Crab are typically observed with 
the standard 5x5 or some other dithering pattern, and this on- 
axis effect is not important for them. For instance, there are only 
19 science windows with Cyg X-1 observed at off-axis angles 
below 0.3°, among all of the 1726 publicly available pointings 
for that source. 

The correctness of the background model provided by the 
background map is crucial for obtaining good PICsIT results. 
As already mentioned before, such a background map cannot 
be prepared with the staring data because the PIE pattern cor- 
responding to the source emission is then included in the back- 
ground pattern and it is not possible to distinguish the signal 
from the background. Tests performed with the maps based on 
the staring data show the Crab and Cen A spectra with a low 
intensity, similar to the noise spectra derived for empty-field ob- 
servations. This does not mean that the staring data must always 
be neglected. Usually only part of the revolution is consumed 
by these observations and the background map can be prepared 
with the rest of the data, corresponding to the dithering mode. 
Even when there is no dithering data for a given revolution, it 
is still possible to prepare background maps using data from 
observations carried out some time before or later In the case 
of the Crab, this solution was tried for Revs. 0039 and 0184. 
Unfortunately, the change in the background pattern after pas- 
sage through the Earth radiation belts appears to be too large 
and the staring-data spectra obtained for these two observations 
were unsatisfactory. The highest quality spectrum for Rev. 0039 
was obtained with the background map of Rev. 0044, showing 
a normal Crab spectral shape but with an about 30% too high 
normalization. The maps from Revs. 0042, 0043, and 0045 pro- 
duced lower quality results, of higher flux and a distorted spec- 
tral shape at energies higher than a 500 keV. This test indicates 
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Fig. 12. Test of the ofF-axis effect. The plots show the ratio be- 
tween the data and a power-law model for Crab observations 
from Revs. 0170-0636 with different selections of off-axis an- 
gles. Model used for each plot was fitted to the spectrum made 
with dithering data. Dithering data comprise all data except on- 
axis data (< 0.3°). For the > 10° selection only the data of Rev. 
0170 were used. 



that the observation closest in time was not necessarily the most 
suitable for preparing a background map. 

Other Crab staring data were collected together with the 
dithering data for the same revolution, allowing us to use the 
default background map. The spectra extracted for both on-axis 
and off-axis staring data show, in general, the same type of be- 
haviour: the flux is too high, and the spectrum is distorted and 
too hard at higher energy. Because the staring observation is 
sometimes performed in different conditions from the dithering 
ones, e.g., at the beginning or end of the revolution, the dither- 
ing background map does not always model sufficiently well 
the background for the entire revolution. This can explain partly 
the observed effect, although an incorrect background map nor- 
mally produces a spectrum that is consistent with the noise level. 
Therefore, the higher flux and harder staring spectra, especially 
those with the source observed on-axis, remain unexplained. 

Figure [13] presents yet another spectral test made for Crab. 
As mentioned in Sect. |2] the wall transparency of the IBIS tele- 
scope allows a quite large fraction of high-energy photons from 
the source to reach the detector when the object is observed at 
large (> 15°) off-axis angles. If this effect is modelled by the PIF 
in its extended version, one can in principle extract spectra for 
these observations. However, the quality of the result depends on 
the actual observing strategy because the source-PIF patterns for 
the science windows used to prepare the background map should 
be sufficiently different from each other to allow for disentangle- 
ment between the source and the background. The far off-axis 
Crab spectrum shown in Fig. [T3] proves that it is possible to ob- 
tain a good result, at least in the energy range up to wSOO keV, 
where no more significantr departure from the normal spectrum 
is seen. This also indicates that the standard response files can be 
used in such a case. Taking into account the fact that for almost 
all objects observed with INTEGRAL, there are more pointings 
with the ofF-axis angle in the 15°^0° range than in the standard 



dithering < 15° range, plenty additional data can be used. This 
technique needs further studies to be carried out together with 
the investigation aimed at finding the most reliable method of 
verifying the applicability of the background map, mentioned at 
the end of Sect. [14] 
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Fig. 13. Spectra extracted for very large off-axis pointings with 
the use of the PIF including the IBIS walls and hopper model. 
Blue crosses show the total PICsIT Crab spectrum obtained with 
the standard PIF for the mask-coded part of the shadowgram. 
Red crosses are the PICsIT Crab spectrum extracted with the 
extended PIF using the data from the Algol observation in Rev. 
0220, when the Crab was 31-39 degrees ofF-axis. 



6.3. Comparison with ttie OS A 7.0 results 

The spectra extracted using the PPDF method have to be com- 
pared with the spectra extracted with the standard PICsIT anal- 
ysis software distributed within the OSA package. Because the 
PIF-based spectral extraction software included in OSA 7.0 is 
still unreliable dFoschini et al.l [20071) . OSA spectra were pre- 
pared from the sky mosaic images. Both methods of creat- 
ing mosaic images were tested: the first applies the OSA tool 
ip_skymosaic, the second utilizes the varmosaic tool distributed 
within the HEASOFT package. The data analysis performed 
using the OSA followed exactly the instructions given in the 
docum ent "PICs IT analysis made easy" prepared by the IBIS 
Team (iFoschinil [2b07a). Mosaic images based on the single- 
event data were created in the 8 default energy bins, cor- 
responding to the channels of the rebinned response matrix 
pics_srmf_grpj)005 .fits. All parameters were set to the default 
values, except for PICSIT JnCorVar, which was set to 1 to take 
into account an uncertainty in the background map. Crab spec- 
tra were extracted for the same sets of science windows as for 
the PPDF method, in order to test the software performance in 
various conditions. 

Figure [14] and Table [5]present a comparison of the OSA and 
PPDF Crab count rates in a broad, low-energy band. For OSA 
this is the second default channel between 252 and 336 keV, i.e. 
the lowest energy channel unaffected by the track events. For 
the PPDF 22-bin spectra, channels 2-5 were merged, giving the 
count rate in the 255-340 keV band. The Crab light curve ex- 
tracted with the PPDF method is remarkably stable and the re- 
sults are much more precise than those obtained with the OSA 
software. The Crab count rates computed by the PPDF method 
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are almost two times higher than those read from the OSA sky 
mosaic images. Since there are so many different elements in 
these two approaches, it is difficult to identify a main reason for 
such a large count-rate deviation. Possible explanation can be: a 
difference between the PIF model and IBIS mask model used 
in the image deconvolution, different background maps, and 
different approach used to calculate the net source-count rate. 
Moreover, the fact that the count rates extracted from mosaic im- 
ages created using ip^skymosaic and varmosaic sometimes dif- 
fer, shows that the mosaic results can be unstable. Although the 
PICSITJnCorVar parameter was set to 1 in the sky image decon- 
volution, the uncertainty in the count rates extracted from the 
mosaic produced with the varmosaic tool are clearly too small 
compared to the scatter of results. 
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Fig. 14. Comparison between the Crab results obtained with the 
PPDF method and with the standard OSA 7.0 software. The up- 
per part shows the Crab light curve in the 255-340 keV band 
extracted with the PPDF method. In the lower part. Crab light 
curves in the 252-336 keV band made with OSA 7.0 are shown. 
OSA mosaic images were obtained with two alternative meth- 
ods, one is the standard OSA tool ip_skyjnosaic and the sec- 
ond is the varmosaic tool distributed in the HEASOFT package. 
Dashed lines show the mean count rate level and dotted lines 
one standard deviation limits. Points showing the ip_sJcy_mosaic 
results are slightly shifted to the right for a better visibility. 



A detailed analysis of the standard OSA 7.0 PlCslT spectra 
is beyond the scope of this paper Nevertheless, we completed 
spectral fitting for about 20 different Crab spectra extracted for 
different observing conditions, which were fitted separately first 
and then together, for both mosaicing tools. In addition, the long- 
exposure spectra from mosaics containing data from many rev- 
olutions were prepared to be compared with PPDF spectra ex- 
tracted for the same data set. The default energy bins defined in 
the OSA software are too wide to allow a detailed spectral analy- 
sis. Therefore, only some basic tests were completed by fitting a 
power-law model to all 8 bins (203-6720 keV) or to the limited 
band with channels 1 and 8 excluded. Spectra of shorter expo- 
sure times, e.g. those corresponding to a single revolution, were 
fitted by retaining the relative normalization as a free parameter. 



Table 5. Parameters of the Crab light curves obtained with the 
PPDF and OSA 7.0 (I - ip_skymosaic, V - varmosaic) methods. 



Method 


Mean 


Std. dev. 


Rel. error 


;r'/NDF 


PPDF 


4.535 


0.155 


3.4% 


12.0/13 


OSA, I 


2.372 


0.422 


17.8% 


65.5/13 


OSA, V 


2.585 


0.473 


18.3% 


15.0/13 



Table 6. Parameters of the Crab spectral fitting obtained for 
dithering spectra extracted with the PPDF and OSA 7.0 (I - 
ip_skymosaic, V - varmosaic) methods. The PPDF spectrum was 
fitted in the 298-1938 keV range, OSA spectra in the 252-1848 
keV range. Flux in the 300-1000 keV band is given in ergs cm"^ 
s Rows 2 and 3 show the OSA results obtained when single 
revolution spectra were fitted together with a free relative nor- 
malization, rows 4 and 5 correspond to the fit made for spectra 
extracted from all-data mosaics. 



Method 


r 


Norm. 


Flux 




Red.^- 


PPDF 


2.07±0.05 


6.2+1.7 


7.8 X 10" 


-y 


0.80 


OSA, I 


2.02±0.16 


2.7±2.5 


4.6 X 10- 


-9 


0.19 


OSA, V 


2.05±0.07 


3.4+1.3 


4.1 X 10" 


-9 


0.97 


OSA, I 


1.98±0.21 


1.6±1.8 


3.3 X 10 


-9 


0.14 


OSA, V 


2.02±0.07 


2.3±0.9 


4.0 X 10- 


-9 


2.39 



In general, the Crab spectral-fitting results were consistent 
for all three methods, i.e. the PPDF and both OSA-based mo- 
saics. The relative normalization of the short-exposure Crab 
spectra was similar for each method, confirming, for example, 
the effect of higher count rates, as measured for on-axis spectra. 
However, these numbers and the spectral indexes are affected by 
significantly more scatter for the OSA results than for the PPDF 
results. The most extreme cases are staring spectra from Rev. 
0102, which have a F value of 3.93, and Rev. 0045 spectra with 
a r of 1.43. Both observations were performed with a relatively 
small Sun off-axis angle and may have been affected by a higher 
and more variable background than usual. For other spectra, F 
values are within the range 1.5-3.1 and the relative normaliza- 
tion parameters are within the range 0.5-1.4. The results of the 
spectral fitting for all dithering spectra together are presented 
in Tabled All spectra fitted in the 300-1000 keV band had a 
slightly smaller spectral index than obtained for the SPI spectra. 
This suggests that the steepness of the PICsIT ARF was over- 
estimated in this energy range. Fluxes obtained with the PPDF 
method are about a factor of two higher than those of the OSA 
spectra, in agreement with the result for a narrower band, pre- 
sented in Table |5] 

PPDF and OSA spectra were also compared for two other 
sources, and the results are shown in Figs. [15] and \TE\ The 
Cyg X-1 observation completed in September 2006 (Rev. 0482) 
was special because the source was then brighter than the Crab. 
Unfortunately, this was a routine Galactic Plane Scan observa- 
tion and the object was always at least 9° off-axis. Therefore, 
the results were not optimal for studies of the physics of this 
Cyg X-1 brightening but, on the other hand, provided a good op- 
portunity to test the OSA performance under special conditions. 
The model plotted in Fig. [15] is a cut-off power-law fitted only 
to the SPI and PICsIT PPDF spectra, in the 100-600 keV band. 
Because this observation took place well after the period cov- 
ered by the data used to prepare the default OSA PICsIT back- 
ground maps, a revised map was prepared for the default OSA 
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energy bins with the data taken from revolution 0482. However, 
both types of PICsIT OSA spectra extracted for this observa- 
tion, i.e. using either the default or new map, were of quite low 
quality. There was no detection (negative count rate) in the first 
and second band of the default-map and new-map spectra, re- 
spectively. The default-map spectrum exhibited large deviations 
from the model, and both of the OSA spectra exhibit high lev- 
els of flux, above that of the PPDF spectrum, whereas the Crab 
results showed the opposite trends. Noise-level studies for an ob- 
servation of similar exposure time (see Sect.jSj indicated reliable 
source detection for energies of up to 600 keV. Thus, the spec- 
tral hardening observed in the PPDF spectrum above that energy 
was only the noise fluctuation. A detailed analysis and discus- 
sion of the INTEGRAL results for the Cyg X-1 flare o bserved in 
September 2006 were presented in lMalzac et alj ( I2008II) . 
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Fig. 15. Comparison between the SPI and PICsIT spectra of Cyg 
X-1 during the September 2006 outburst. The model was fitted 
to the SPI and PICsIT PPDF spectra in the 100-600 keV band. 



The third object, for which the OSA PICsIT spectra were 
tested was Cen A, to check whether is was possible to study a 
source much weaker than Crab or Cyg X-1. The same data as 
in the case of the calibration studies were used, and the OSA 
mosaic was made with the varmosaic tool only. The first two 
channels of the OSA PICsIT spectrum are in agreement with 
the model, assuming a relative normalization similar to the Crab 
case. Nevertheless, the OSA spectrum is much less precise than 
that extracted with the PPDF method and provides no additional 
information to that derived using the SPI spectrum. 

6.4. PICsIT detections 

We review PICsIT spectra extracted for objects of various types. 
Although most of these results are only preliminary and were 
obtained without a careful analysis and data selection, they illus- 
trate the type of information that PICsIT can provide. The num- 
ber of objects with relatively strong emission within the PICsIT 
energy range between several hundreds of ke V and several Me V 
is very small. The dominant class of sources are Galactic binary 
systems with a black hole candidate. Three other sources de- 
tected until now were the Crab and PSR B1509-58 pulsars and 
a single AGN, Cen A. Several objects were detected only dur- 
ing outbursts of several days or weeks because their persistent 
emission was too weak for PICsIT. Currently, there are 12 firmly 
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Fig. 16. Comparison between the SPI and PICsIT spectra of Cen 
A observations in 2003-2004. The model was fitted to the SPI 
and PICsIT PPDF spectra. 



detected sources (with > 6-sigma detections in the 277-361 keV 
band) among about 30 tested. 

Figure[T7]illustrates the scale of the problem for a detection 
of emission in the soft y-ray domain. Cosmic-ray induced inter- 
nal background of the detector is usually orders of magnitude 
higher than the radiation of the brightest sources. In the case of 
PICsIT, the background is a factor of 100 stronger than the Crab 
emission around 300 keV, and 1000 times stronger at several 
MeV. Due to the coded-mask technique and dithering observa- 
tion strategy, it remains possible to model the background to a 
precision allowing objects many times weaker than Crab to be 
detected, provided that there is sufficient exposure time. Crab 
was observed only during one or two revolutions twice a year, 
and the total exposure time of about 1 Ms provides a clear de- 
tection up to 2 MeV. Cyg X-1 was observed far more frequently 
but its steeper spectrum allowed us to detect the source only up 
to about 700 keV. Because more data will become public from 
the INTEGRAL Key Programme observations of the Cygnus re- 
gion, future analysis should yield a spectrum reaching at least 1 
MeV, according to the noise limits shown in Fig. [17] 
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Fig. 17. PICsIT spectra of the two brightest objects. Crab and 
Cyg X-1, compared to the PICsIT background spectrum and 3- 
sigma upper noise limits for 1, 2 and 14 Ms of eff'ective exposure 
time. 
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The spectra of three other persistent sources that can be well 
studied with PICsIT are presented in Fig. [18] Due to a long ex- 
posure time, the results obtained for the microquasar GRS 1758- 
258 were of overhelmingly higher quality than any obtained for 
this object with any other, past, or contemporaneous, instrument. 
As data continues to arrive in the next few years for this Galactic 
Bulge source it will be possible to study the emission proper- 
ties up to about 600 keV, with the possibility of subdividing the 
data according to the object state. The other bright microquasar, 
GRS 1915-H105, would need a more sophisticated analysis, but 
preliminary results indicate rather unambiguously a harder spec- 
trum in the soft y-ray band than for GRS 1758-258, and a cut-off 
at around 500 keV. The third object with the spectrum shown in 
Fig.[ll IE 1740.7-2942, exhibits weaker emission in the PICsIT 
energy band, and more data and a more careful analysis will 
be needed to obtain a high-quality spectrum. Nevertheless, the 
final spectrum should not differ significantly from the current 
one, showing yet another emission pattern with a hard spectral 
slope. The other spectrum of Fig. [18] was extracted for the posi- 
tion of the Galactic Centre radio-source Sgr A*. Taking into ac- 
count that no strong hard X-ray emission originates in this point 
(ISGRl observations of IGR J 17456-2901 reveal a 7 mCrab 
source in the 60-150 keV band (iKuulkers et al.ll2007l) ). it is un- 
surprising that the PlCslT spectrum is fully consistent with the 
noise spectrum. The expected exposure time of in excess of 20 
Ms at the end of the INTEGRAL mission will allow us to derive 
the most reliable ever upper limit to the emission from the centre 
of our Galaxy. We emphasize that for both IE 1740.7-1-2942 and 
Sgr A* located in a region of str ong diffuse emission observed 
in the 200-500 keV range by SPI (iBouchet et al ."2008). PICsIT, 
due to its high angular resolution, is the only instrument capable 
of resolving point sources without the ambiguity affecting SPI 
and Suzaku/HXD results. 

Besides the bright persistent soft y-ray emitters mentioned 
above, several objects are too weak to be detected by PICsIT in 
the normal state but bright enough during relatively long periods 
of outbursts. Five sources of this type have been detected with 
PICsIT. Figure[T8]presents the spectrum of one of these sources, 
XTE J1550-564. The other four are black hole candidates: GRO 
J1655-40, GX 339-4, IGR J17497-2821, and SWIFT J1753.5- 
0127. The last spectrum shown in Fig. [18] illustrates detection 
of emission from the PSR B1509-58 pulsar system. This source 
was d etected by PICsIT with the standard OS A software dKuipeij 
120051) . although no details were then given. 



6.5. Future prospects 

The spectral extraction method presented in this paper provides 
valuable and reliable results, as shown in the previous subsec- 
tions. Notwithstanding, possible ways of improving the PICsIT 
spectra remain, which would increase the sensitivity of the in- 
strument and the stability of the results. The first possibility is 
application of a spectral extraction also to the multiple events. 
Although several times fewer photons are registered in multiple 
than single events below 1 MeV, the instrumental background 
for multiple events is also several times lower. This can allow 
the sensitivity to be improved by about 20% above 300 keV. 
More careful modelling of the background and more careful se- 
lection of the data will produce more stable source count-rates, 
smoother spectra, and more precise light curves. This goal can 
be achieved by developing tools to verify the adequacy and sta- 
bility of the background during a single revolution or a sequence 
of observations. Currently, a simple way of achieving a more ac- 
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Fig. 18. Examples of PICsIT spectra extracted for several 
Galactic objects. Spectrum extracted for the position of Sgr A* 
presents situation, when there is no detection. Spectra in four 
upper panels are compared with the 3-sigma upper noise lim- 
its derived for an effective exposure time of 10 Ms. Spectrum 
of XTE J 1550-564 measured during outburst is compared with 
the 3-sigma upper noise limits computed for an effective expo- 
sure time of 200 ks and the PSR B1509-58 spectrum with the 
3-sigma upper noise limits for an effective exposure time of 2 
Ms. 



curate background model is to prepare two or three background 
maps for a given revolution. 

Further improvements in sensitivity, of up to a factor 2, can 
be achieved with the development of an efficient tool to verify 
the background-map applicability to the data of other revolu- 
tions. This would allow us to extract the results for sources ob- 
served in staring mode or at very large off-axis angles, when the 
standard map produced for a given revolution does not allow us 
to distinguish between the source and the background. In this 
way, the total exposure time for many objects could be increased 
by a factor of 2^, which would be important for extragalatic 
objects observed rarely with INTEGRAL. 



7. Summary and conclusions 

We have developed a comprehensive method for extracting 
PICsIT spectra from spectral-imaging single-event data, which 
enables us to derive information to the physical limitations of the 
instrument. This was achieved mainly by application of a correct 
description of the Poisson-distributed data and careful modelling 
of the detector background. 

We adopted a technique that handles the Poisson probabil- 
ity density functions in a Bayesian manner for the first time in 
extracting the count rates from a y-ray detector The effective- 
ness of this technique in the 'low number of counts' regime and 
its superiority over other methods was proven with exhaustive 
tests. 
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In our spectral extraction technique we employed a back- 
ground model that was as close as possible to the data because 
it consisted of a map of the mean count rate measured for each 
pixel during a given object observation. This approach is justi- 
fied for a coded-mask instrument operating in a dithering mode, 
especially when the background is orders of magnitude stronger 
than the source emission. 

The detection reliability was verified by detailed studies of 
fake source spectra extracted from empty field observations. 
This approach provided the most reliable estimate of the noise 
level possible, which accounted for true background fluctuations 
and other systematic effects related to the instrument and spec- 
tral extraction method. We recommend similar tests for other 
instruments, because the standard analytical calculation of the 
detection limit can often underestimate the true noise limit. 

We tested the performance of the new PICsIT spectral ex- 
traction method by comparing our results with those from other 
INTEGRAL instruments. Using the standard PICsIT response 
files, we found good agreement between the spectral fitting re- 
sults for data from PICsIT and both SPI and ISGRI, obtained 
for the Crab, Cen A, and Cyg X-1. PICsIT appears to be well 
calibrated, although some tuning of the response files may be 
needed below 300 keV and above 1 MeV. 

Confronted with the standard OSA software results, count 
rates computed with the new method are far more stable and 
about a factor of two higher, in close agreement with results de- 
rived for data from other INTEGRAL instruments. Crab spectra 
extracted with both OSA and the new method are described well 
by a model of the same spectral slope. However, for weaker ob- 
jects and for observations made at large off-axis angles, we were 
unable to obtain high quality spectra with the standard software. 

The new spectral extraction method was applied to date for 
about 30 objects observed with INTEGRAL. Eight new sources 
were added to the four already reported, which were detected 
with the standard PICsIT software. A refined analysis and an ex- 
panding data set should provide at least several detections more. 
Due to a very long exposure time, the spectrum obtained for the 
microquasar GRS 1758-258 reached an unprecedented quality 
and could be extended to higher energies up to about 700 keV. 
Preliminary spectra extracted for several other sources, in partic- 
ular those observed during outbursts, also provide information 
unachievable for any other instrument currently operating. 

Despite its unique advantages, PICsIT has a limited sensitiv- 
ity when compared to the OSSE or BATSE detectors. The main 
reason for this is the smaller thickness of the detector layer, 3 
cm for PICsIT, compared to 17.8 cm for OSSE and 7.6 cm for 
BATSE. In addition, the PICsIT ability to detect photons is de- 
creased by the pixellated structure of the detector and the logic 
discriminating single and multiple events. Nevertheless, it seems 
that a coded-mask instrument is the most appropriate solution for 
soft y-ray astronomy, because this technique offers at the same 
time a large field of view, good angular resolution and direct 
measurement of the background when dithering observations are 
used. 
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Appendix A: Tests of spectral extraction methods 
applied to the ISGRI detector 

The main aim of the tests presented here was an examination 
of the behaviour of different spectral extraction techniques in 
a situation when the source emission is very weak. The stan- 
dard OSA software for ISGRI spectral extraction is the PIF- 
based iLspectra_extract tool, using th e maximum likelihoo d 
technique with Poissonian distribution dGoldwurm et al.ll2003h . 
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A detailed description of this technique will be published soon 
(A. Goldwurm and A. Gros, in preparation). Besides the OSA 
7.0 standard method and the PPDF method appUed to the ISGRI 
spectral extraction, five alternative methods were checked. Two 
of them also represent a Bayesian approach with the direct han- 
dling of probability density functions, where the Poisson dis- 
tribution is replaced by the Gaussian PDF (hereafter GPDF). 
The first one holds the Gaussian width parameter cr equal to the 
square root of the number of counts given by the model, whereas 
the second uses instead a cr equalling the square root of the sim- 
ulated number of counts. As in the case of the PPDF method, 
the final outcome of these two techniques is the probability den- 
sity distribution for the source net count rate. Then, two standard 
statistics methods were controlled, where the test statistic was 
either a (Gaussian likelihood) or C statistic (Poisson likeli- 
hood). For these two techniques the result collected for a single 
science window consisted of only two numbers: the source count 
rate minimizing the (or C) and its uncertainty corresponding 
to a change of the statistic by a certain value. In the case of the;\f^ 
statistic, both variants of the Gaussian width parameter (from the 
model or from data) were checked. When the number of counts 
simulated for a given pixel was zero, the 'data variant' cr was set 
to be 1. 

Final results for the methods producing probability density 
distributions were computed using the product of all individual 
PDFs obtained for single science windows constituting a given 
data set (simulation run). OSA spectra were summed with a ded- 
icated tool spe_pick (included in the OSA software). In the case 
of methods based on x^ and C statistic, the results for a single 
science window were computed over a grid of the source count 
rate and background normalization values. This was more time 
demanding, but allowed us to avoid potential problems with the 
convergence of the fitting procedure when a global minimum 
(maximum) is searched for. Uncertainties in a single science 
window result were determined by finding the width of the confi- 
dence region projected onto the source count-rate axis, for which 
the x^ or C statistic value changed by 1.0, i.e. 1-signi a confi- 
dence level for one interesting parameter dCashl [19791) . Then, 
the count rates for larger data sets (e.g. one or all simulation 
runs) were calculated as a mean weighted by these uncertainties. 
This approach is, of course, not justified for the C statistic be- 
cause the likelihood obtained for a single science window is still 
not well approximated by the Gaussian distribution. In conse- 
quence, the standard weighted mean is inadequate for an asym- 
metric Poisson distribution characterizing 'low count rate' data, 
where the position of the density maximum is usually far from 
the position of the mean. Because the maximum likelihood is of- 
ten peaked at zero and the lower and upper limits are located at 
quite different distances from the maximum, the standard devia- 
tion error used for computing the weighted mean was set to the 
arithmetic mean of the lower and upper errors. 

Simulations were performed for the ISGRI detector with 
the aim of comparing the PPDF method with the standard 
OSA ISGRI results and because the method presented here was 
planned to be used with the high-energy ISGRI data. The obser- 
vations of the Cassiopeia field made in Revs. 0384-0395 served 
as a simulation template. They included 428 science windows 
with about 880 ks total exposure time. A fake source of 1 mCrab 
strength was assumed to be at the Cas A position. The standard 
OSA 7.0 PIFs computed for Cas A and background maps were 
used to simulate the shadowgrams in thirteen energy bands listed 
in Table lA.fl The simulated source and background count rates 
are also presented, together with the signal-to-noise ratios (S/N) 
for 1 Ms exposure time. The noise was estimated simply to be 



Table A.l. Initial conditions of tests done with the simulated 
ISGRI data. The signal to noise ratio, S/N, was calculated for a 
1 Ms effective exposure time. 



Energy range Source count rate Bkg. count rate S/N 



[keV] 


[cps] 


[cps] 




14.0_17.8 


0.035 


31.1 


6.3 


17.8-21.6 


0.037 


35.8 


6.2 


21.6-25.4 


0.032 


39.7 


5.1 


25.4-31.2 


0.035 


44.1 


5.3 


31.2-37.9 


0.026 


35.7 


44 


37.9-45.6 


0.022 


34.9 


3.7 


45.6-55.1 


0.021 


48.9 


3.0 


55.1-68.5 


0.021 


92.0 


2.2 


68.5-85.8 


0.013 


68.5 


1.6 


85.8-112.6 


0.013 


65.1 


1.6 


112.6-166.2 


0.0083 


85.6 


0.9 


166.2-276.3 


0.0023 


105.2 


0.22 


276.3^62.5 


0.00033 


69.9 


0.04 



the ratio of the square root of background counts to the exposure 
time. 

Because the simulated count rates were too low to achieve a 
3-sigma detection in the higher energy bands, simulations were 
repeated 20 times, resulting in 17.6 Ms of total simulated expo- 
sure time. In addition, to check the case with conditions more 
relevant to PICsIT, another set of simulations with the back- 
ground count rates assumed to be 10 times higher than those 
of ISGRI was performed. 

Figure [ATl presents the main result of the tests. For the stan- 
dard background level, both PPDF and OSA 7.0 methods repro- 
duce correctly the source count rates with at least a 3-sigma sig- 
nificance for energy bands below 166 keV. The PPDF method 
performs more reliably, producing results very close to the as- 
sumed values in all bands below 40 keV. Above that energy, 
PPDF count rates are, on average, closer to the true value than 
the OSA 7.0 count rates. In the case of a test made with a back- 
ground level higher by a factor of 10, the differences between 
the PPDF and OSA 7.0 results are small in the 18-55 keV en- 
ergy range, although, at higher energies the PPDF method ap- 
pears to estimate the count rate more reliably (note that there is 
still insufficient data to obtain a detection at higher energy). 

Since alternative methods were not expected to produce good 
results at low count rates, their behaviour was tested first for a 
wider range of simulated count rates. Figure |AT2l shows the re- 
sults of this test for the second energy band (17. 8-2 1.6 keV) with 
the count rates decreasing from 10 to 0.02 cps. The GPDF and 

statistic methods using the model cr parameter diverge from 
the true result by a factor of 2 or more starting from 0.2 cps, pro- 
ducing too high count rates. When the data cr parameter is used 
instead in these two approaches, the reproduced count rates are 
too small by a factor of > 2 already at 0.5 cps. The C statistic per- 
forms much better than the methods based on the Gaussian dis- 
tribution: deviations from the true value are observed below 0.5 
cps but they remain small (< 20%), except for the lowest tested 
count rate (0.02 cps). Below 0.05 cps, the count rates reproduced 
by the GPDF, the x^ statistic (both with the model cr) and the C 
statistic all saturate at different levels. We checked with an ad- 
ditional test that these levels are proportional to the background 
level: for a 10 times higher background they are about 3 times 
higher, in accordance with the expected increase in the standard 
deviation, equal to the square root of the total count number. 
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Fig.A.l. Results of tests (all 20 simulation rans merged to- 
gether) for the two best spectral extraction methods: the new 
PPDF technique and the standard ISGRI spectral extraction soft- 
ware released in OSA 7.0. Upper panel: ratios between the ex- 
tracted and assumed count rates for the normal ISGRI back- 
ground level. Lower panel: ratios obtained for a background 
10 times higher. The energy bands of OSA results are shifted 
sUghtly to the left for a better visibility. 



Figure IA3l presents again the comparison between the PPDF 
and other methods, this time for all tested energy bands and 
for the count rates corresponding to the 1 mCrab source ob- 
served with ISGRI. At this count rate level, only the PPDF tech- 
nique provides the correct results. The deviations observed for 
the other methods are correlated with the increasing background 
level (see Table |A. ft . At low energy (below 20 keV), when the 
number of counts simulated for a given pixel is often equal to 
zero (due to the low energy thresholds set for the pixels), both 
methods using the Gaussian approximation (GPDF and x^) fail 
completely. 
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Fig.A.l. Comparison between the results of the best method 
(PPDF) and the results obtained with the use of several alter- 
native methods. Tests were done for the 17.8-21.6 keV energy 
band and the standard ISGRI background. One simulation run 
(428 science windows) was done for the count rates above 0.2 
cps. For 0.2, 0.1, 0.05 and 0.02 cps there were 2, 5, 10 and 10 
simulation runs merged together, respectively. 
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Fig.A.3. Comparison between the results of the best method 
(PPDF) and the results obtained with the use of several alter- 
native methods. Tests were done for the standard ISGRI back- 
ground. 



Figure |A3I compares the results of the PPDF and C statistic 
methods applied to single science windows and merged later in 
the way described above, and the results obtained when all the 
data from the entire simulation run (428 science windows) were 
extracted at once with these methods. Merging the data is pos- 
sible because simulations are always performed with the same 
background level assumed and all 428 data sets can be treated 
as one set, corresponding to the detector with 428 times more 
pixels. Contrary to the standard procedure, the background nor- 
malization was found once for all science windows. For both 
higher and lower S/N bands (low and high energy, respectively), 
both PPDF approaches give similar results, completely consis- 
tent, on average, with the assumed count rate. Direct applica- 
tion of the PPDF count rate extraction to the data of all science 
windows at once leads to slightly larger final errors, because 
of a broader two-dimensional PDF associated with a broader 
background count rate distribution in a larger simulation set. It 
seems more reliable to integrate the two-dimensional PDF over 
the background normalization factor for each science window 
separately because it allows a more careful control of back- 
ground fluctuations. Since in a real situation the instrumental 
background varies quickly, this strategy should be followed and 
one has to work on possibly short-time data sets before merging 
the results from a longer observation period. 

On the other hand, the C statistic technique produces the cor- 
rect results only when applied to a larger data set. In this case, the 
maximum likelihood method of course finds a correct position of 
the PDF maximum, as shown in Fig. I A. 41 where the PPDF max- 
ima are marked for comparison. Nevertheless, it appears to be 
rather inefficient for low S/N data, as shown in the lowest panel 
of Fig. [A4l This is a consequence of using the maximum, which 
is quite often located far from the gravity centre of the asymmet- 
ric Poisson PDF. In principle, when the count rates are not too 
small, the C statistic applied to all data at once can be used, pro- 
vided that there is no need to update the result when a new data 
set arrives. For such an update, the PDF methods storing all dis- 
tributions derived for single science windows are faster because 
in this case one has only to extract the PDFs for the new data and 
then merge them with the older results. Also any change in the 
data selection, e.g. for the source variability studies, makes the 
method based on C statistic rather impractical. 
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Fig. A.4. Results of two variants of the C statistic method com- 
pared to the corresponding resuhs for the PPDF method. The 
comparison is made for two energy bands representing high (at 
low energy) and low (at high energy) S/N ratios. In each panel 
the points shifted to the right illustrate the results of a given 
method applied to each science window separately and then 
merged by multiplying all PDFs (PPDF) or by computing the 
weighted mean (C statistic). Points shifted to the left represent 
the variant of applying either method to the total data set at once. 
Open diamonds indicate the position of the maximum of the 
PPDF distribution. Dashed lines show the count rate assumed 
in simulations. 



Besides the performance tests described already, the two 
most reliable methods were in addition checked with respect to 
the scatter in the results and the correctness of the derived un- 
certainties. These tests were done using the results of 20 simu- 
lation runs and merging them in a random way to obtain a new 
simulation run, again with 428 science windows. For each new 
run, each science window was selected randomly from one of 
20 runs, to preserve the original science-window content of the 
simulation set. This procedure was repeated 500 times, to cal- 
culate arithmetic mean of the standard deviations &, and the 
fraction of results containing the true count-rate value within the 
uncertainty limits. The complete results of these tests are listed 
in Table |A.2l for the first 1 1 energy bins and for three energy bins 
(4, 7, and 10), for the default and 10 times higher background, 
respectively. 

Both methods were examined with the so-called 'frequen- 
tist coverage' test, where the fraction of results covering the 
true count rate value by error bars was computed. In case of the 



Table A.2. Uncertainty tests done for the OSA 7.0 and PPDF 
spectral extraction methods, is summed over all 500 shots, 
& is an arithmetic mean of the standard deviation errors, /lo- is 
a fraction of shots covering the true count rate value inside the 
error limits and /, is a similar fraction of shots computed for 
the shortest credible interval errors for the PPDF method. Three 
rows at the bottom present the results of the same test made for 
the background assumed to be 10 times higher than the default 
one. 



Energy range 




OSA 7.0 






PPDF 




[keV] 




(T 


ficr 




cr 




fs 


14.0-17.8 


All 


0.0114 


0.69 


425 


0.0089 


0.73 


0.72 


17.8-21.6 


566 


0.0087 


0.65 


537 


0.0084 


0.66 


0.64 


21.6-254 


467 


0.0086 


0.67 


443 


0.0084 


0.71 


0.71 


25.4-31.2 


478 


0.0090 


0.68 


503 


0.0087 


0.68 


0.68 


31.2-37.9 


481 


0.0081 


0.72 


459 


0.0078 


0.73 


0.71 


37.9-45.6 


513 


0.0080 


0.67 


569 


0.0076 


0.65 


0.63 


45.6-55.1 


489 


0.0094 


0.70 


486 


0.0085 


0.72 


0.68 


55.1-68.5 


508 


0.0130 


0.69 


443 


0.0105 


0.74 


0.56 


68.5-85.8 


447 


0.0112 


0.71 


344 


0.0088 


0.78 


0.65 


85.8-112.6 


462 


0.0109 


0.71 


285 


0.0083 


0.82 


0.61 


112.6-166.2 


460 


0.0125 


0.70 


459 


0.0091 


0.71 


0.75 


25.4-31.2 


537 


0.0284 


0.71 


348 


0.0225 


0.78 


0.66 


45.6-55.1 


487 


0.0299 


0.67 


435 


0.0221 


0.72 


0.75 


85.8-112.6 


513 


0.0344 


0.65 


367 


0.0207 


0.76 


0.88 



PPDF method, there are many ways to determine the uncertainty, 
so-called 'credible intervals' in Bayesian nomenclature, because 
the PDF associated with the result is asymmetric. Here, three 
variants were tested. The first one was a standard deviation de- 
fined as the square root of the variance determined by integrating 
is - sf- (where s is the source count rate) weighted by the PDF, 
i.e. the second central moment of the PDF. The second variant 
was the shortest interval containing the desired fraction of the 
PDF integral, 68.3% for l-cr error. The third option was the cen- 
tral credible interval, found as limits above and below which the 
PDF integrals were equal to (l-0.683)/2. 

The standard deviation and central credible interval ap- 
proaches gave almost identical results and thus only the former 
is shown in Table IA.2I They tend to overestimate slightly the 
uncertainty, whereas the shortest interval approach seems to un- 
derestimate the uncertainty. Nevertheless, all of these variants 
produce errors that are, in general, in agreement with the 0.683 
coverage probability taking into account about 5% deviation ex- 
pected for 500 simulations performed. The OSA 7.0 uncertain- 
ties appear to pass the coverage test more successfully than those 
obtained with the PPDF technique. However, this happens at the 
cost of much larger error values (measured by a-) and a larger 
scatter in the results. Plainly the PPDF method is more precise 
than the method implemented in the OSA software. 

Yet another comparison between the results of OSA 7.0 and 
the PPDF methods is shown in Fig. IA.5I Progress in the deter- 
mination of the source count rate is presented for three energy 
bands. The PPDF produces overestimated count rates when there 
is insufficient information in the data, as a consequence of broad 
but non-negative PDFs. On the other hand, the OSA 7.0 soft- 
ware has some tendency to underestimate count rates when the 
data set remains too small for a detection. With the increasing 
amount of data both approaches exhibit normal random oscil- 
lations around the true value, with the amplitude of variations 
decreasing with the data set size. As already demonstrated in 
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Table IA.2I the scale of oscillations is clearly smaller in the case 
of the PPDF method, converging quite rapidly to the true value. 
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Fig.A.5. Results of tests for the two best spectral extraction 
methods, PPDF and OSA 7.0. Progress in the accuracy of the 
extracted count rates with the increasing number of simulated 
science windows is presented for several energy bands. 



Extensive tests performed with simulated data showed that 
only two count-rate extraction methods work correctly for data 
of very low signal-to-noise ratio. The method handling the 
Poisson probability functions can be treated as a reference, be- 
cause it allows us to obtain the most precise and quickly converg- 
ing results. The only limitation of this method is a long compu- 
tation time when many parameters have to be treated at once. 
However, in the case of high energy astrophysics the number 
of sources in the field of view is always small and this limita- 
tion poses no serious problem. The technique implemented in 
the OSA software also produces very good results and can be 
recommended for general use whenever computation speed and 
robustness is needed. Among the other methods tested, only that 
based on the C statistic performs relatively well when the signal 
is not extremely weak. Because it can produce systematically 
underestimated results, its performance in a given application 
should be tested with simulated data. 
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